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Abstract: Most remote sensing studies assess the desertification using vegetation monitoring method. But 
it has the insufficient precision of vegetation monitoring for the limited vegetation cover of the 
desertification region. Therefore, it offers an alternative approach for the desertification research to assess 
sand dune and sandy land change using remote sensing in the desertification region. In this study, the indices 
derived from the well-known tasseled cap transformation (TCT), tasseled cap angle (TCA), disturbance 
index (DI), process indicator (PI), and topsoil grain size index (TGSI) were integrated to monitor and assess 
the desertification at the thirteen study sites including sand dunes and sandy lands distributed in the 
Mongolian Plateau (MP) from 2000 to 2015. A decision tree was used to classify the desertification on a 
regional scale. The average overall accuracy of 2000, 2005, 2010 and 2015 desertification classification was 
higher than 90%. Results from this study indicated that integration of the advantages of TCA, DI and TGSI 
could better assess the desertification. During the last 16 years, Badain Jaran Desert, Tengger Desert, and 
Ulan Buh Desert showed a relative stabilization. Otindag Sandy Land and the deserts of Khar Nuur, Ereen 
Nuur, Tsagan Nuur, Khongoryn Els, Hobq, and Mu Us showed a slow increasing of desertification, whereas 
Bayan Gobi, Horgin and Hulun Buir sandy lands showed a slow decreasing of desertification. Compared 
with the other 11 sites, the fine sand dunes occupied the majority of the Tengger Desert, and the coarse 
sandy land occupied the majority of the Horgin Sandy Land. Our findings on a three or four years’ 
periodical fluctuated changes in the desertification may possibly reflect changing precipitation and soil 
moisture in the MP. Further work to link the TCA, DI, TGSI, and PI values with the desertification 
characteristics is recommended to set the thresholds and improve the assessment accuracy with field 
investigation. 
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1 Introduction 


Desertification poses serious ecological, environmental, and socio-economic threats to the world. 
It directly affects 2.5x10* people and one-third of the Earth’s surface (over 4.0x107 km?; UNCCD, 
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2016a). According to the definition put forward by the United Nations Convention to Combat 
Desertification (UNCCD) in 1994, desertification is land degradation in arid, semi-arid, and dry 
sub-humid areas due to climatic variation and human intervention. The Asian continent is occupied 
by approximately 1.7x10’ km? of arid, semi-arid, and dry sub-humid lands and desertification 
manifests itself in many different forms, such as the expanding deserts in China, Mongolia, India, 
Iran, and Pakistan (UNCCD, 2016b). As the second largest plateau in Asia, the Mongolian Plateau 
(MP) consists primarily of the Inner Mongolia Autonomous Region of China (IMAR, a total area 
of 1.18x10° km?) and the entire territory of Mongolia (a total area of 1.57x10° km’; Fang et al., 
2015). Approximately 90% of Mongolia’s territory is occupied by hyper-arid, arid, semi-arid, and 
dry sub-humid areas, and roughly 72% of the country is vulnerable to desertification (Yu et al., 
2013; Eckert et al., 2015). In the IMAR, 60% of the land area is occupied by arid, semi-arid, and 
dry sub-humid areas and suffers from desertification (Ci and Wu, 1997). Over the past few decades, 
rapid growth of the population and livestock, urbanization, and mining have exacerbated the 
difficulty in mitigating desertification faced by the government of China and Mongolia. Therefore, 
there is an urgent need to understand and control the desertification in terms of its status, trend, and 
drivers. The most effective way to achieve this is to monitor and map the desertification processes 
(Xiao et al., 2006; Albalawi and Kumar, 2013). 

Traditional desertification monitoring approaches from a site-specific perspective make it 
difficult to understand the spatial extent of the desertification (Collado et al., 2002; Lam et al., 
2010). Remote sensing, with the ability to rapidly and reliably collect data over wide areas, has 
been applied to monitor and assess desertification and its dynamic processes (Cui et al., 2006; 
Albalawi and Kumar, 2013). Changing vegetation and land use are the most common indicators 
used by many scholars to monitor desertification (Albalawi and Kumar, 2013). Generally, changes 
in vegetation are the most direct indicator of desertification (land degradation), and can reflect the 
dynamic processes of desertification (Yang et al., 2005; Xu et al., 2009). Image-based vegetation 
indices such as the normalized difference vegetation index (NDVDJ), soil adjusted vegetation index 
(SAVI), modified SAVI (MSAVI), enhanced vegetation index (EVI), vegetation temperature 
condition index (VTCD), and optimized vegetation index (OPVI), were most commonly used for 
monitoring the vegetation conditions in desertification assessment (Cui et al., 2006; Xu et al., 2009; 
Yu et al., 2013). In addition, some indices that can reflect spatial heterogeneity, such as the land 
surface temperature (LST), land surface albedo, land surface water index (LSWI), moving standard 
deviation index (MSDIJ), and rain use efficiency (RUE), are also used to assess desertification. In 
Mongolia, trend analysis of MODIS NDVI time series data is used to detect areas of vegetation 
change to monitor land degradation and regeneration areas (Eckert et al., 2015). The combination 
of NDVI and field-surveys provides an effective way to monitor desertification on the Mongolian 
steppes (Sternberg et al., 2011). MODIS NDVI cross-checked with the meteorologically-derived 
aridity index from 2000 to 2012 identifies a recent contraction of the Gobi in East Asia (Sternberg 
et al., 2015). Compared with NDVI and EVI, MSAVI2 derived from multi-temporal SPOT-4 
VEGETATION data from 1998 to 2001 is found to estimate vegetation biomass and monitor 
vegetation degradation well on the Mongolian desert steppes (Javzandulam et al., 2005). The VTCI 
from multi-temporal Landsat-5 images has been used to predict the distribution of vegetation and 
efficiently monitor desertification in Bulgan, Mongolia (Yu et al., 2013). The NDVI, MSDI, and 
land surface albedo from Landsat images have been used in conjunction to quantitatively assess the 
desertification of the Ordos Plateau, China in 1980, 1990, and 2000 (Xu et al., 2009). The topsoil 
grain size index (TGSI) is proposed based on the field spectral reflectance measurements and 
laboratory physical analysis of topsoil grain composition to detect desertification in Siziwang 
Banner, IMAR (Xiao et al., 2006). Integrated desertification assessment conducted by combining 
NDVI, TGSI, and land surface albedo from Landsat images indicated a difference of 87% in the 
desertification area in Hogno Khaan nature reserve in central Mongolia between 1990 and 2011 
(Lamchin et al., 2016). The total accuracy of all Landsat MSS-derived vegetation indices (NDVI 
and SAVI) for desertification mapping of the Horgin Sandy Land, IMAR, is approximately 48% 
(Bremborg, 1996). The MSAVI variation is not sufficient to assess the land degradation process, 
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which indicates that the vegetation greenness alone is not a major indicator of land degradation 
assessment in the Gobi-Altai province in western Mongolia (Vova et al., 2015). Land use and land 
cover change (LUCC) has been accepted widely as a primary driver for ecosystem degradation 
worldwide and become a key research topic in desertification. A full-year 2000 time-series of SPOT 
VEGETATION images with 1 km spatial resolution has been used to produce a land cover map to 
detect sparse vegetation as an indicator of areas at risk of desertification in North China, including 
the Gobi desert (Huang and Siegert, 2006). Based on a GIS approach, researchers performed land 
desertification monitoring and assessment in Yulin, Northwest China using spatial patterns from 
land cover classes, Landsat NDVI, and the prevailing wind direction (Zhang et al., 2008). The 
status, rate, and causes of the development of desertification in the Mu Us Sandy Land, northern 
China, has been analyzed based on land cover change from the 1950s to 1990s (Wu and Ci, 2002). 
The spatial and temporal dynamics of land use and land cover associated with land degradation and 
desertification in the Mu Us Sandy Land, China have been analyzed using NDVI and albedo 
through the change vector analysis technique (Karnieli et al., 2014). Monitoring changes in 
vegetation and land use is not the only method available to quantify desertification (Kawamura and 
Akiyama, 2010; Albalawi and Kumar, 2013). Some scholars criticize using vegetation cover to 
detect desertification (Shafie et al., 2012). It is also noteworthy that in the Gobi desert, the NDVI 
must be verified with a ground-based investigation due to the insufficient precision of the NDVI 
value for such limited vegetation cover (Lam et al., 2010; Sternberg et al., 2011; Zhao et al., 2015). 
Currently, it is not feasible to obtain field data and carry out ground-based investigations due to the 
remoteness and very large area of the MP (Sternberg, 2012). 

Sandy desertification is one of the main forms of land degradation in the MP, especially in 
northern China (Elhadi et al., 2009; UNCCD, 2016b). Desert area covers 14.3% of Mongolia (Yang 
et al., 2004), and 40.0% of the IMAR in China (John et al., 2008); this has been steadily expanding 
since the 1950s (Eltahir et al., 2009). The inclusion of sand dune encroachment in desertification 
monitoring would greatly improve the accuracy in the prediction criteria of risk-prone areas (Lam 
et al., 2010; El-Magd et al., 2013). It has been found that sand dune encroachment can be effectively 
monitored through comparison of multi-temporal satellite images based on visual interpretation; 
differences in spectral brightness, land surface albedo, and temperature; image subtraction; and 
supervised and unsupervised classification (Hu et al., 2002; Yao et al., 2007; Hereher, 2010; Lam 
et al., 2010; Hermas et al., 2012; Hugenholtz et al., 2012). However, compared with several discrete 
spectral bands or vegetation indices such as NDVI, MSAVI and normalized difference moisture 
index (NDMI), the tasseled cap brightness (TCB), tasseled cap greenness (TCG) and tasseled cap 
wetness (TCW) derived from the tasseled cap transformation (TCT), which synthesizes the 
information contained in several discrete spectral bands into a few components and can associate 
remotely sensed spectral features with the physical characteristics of land use and land cover classes 
can yield slightly better results (Jin and Sader, 2005; Lozano et al., 2007; Powell et al., 2010; 
Czerwinski et al., 2014; Liu et al., 2016). 

Recently, the different combinations of the tasseled cap indices such as the disturbance index 
(DI) and tasseled cap angle (TCA) have been used in forest disturbance detection and land cover 
change analysis (Masek et al., 2008; Powell et al., 2010; Gómez et al., 2011; Baumann et al., 2014). 
The DI is a linear combination of the normalized TCB, TCG, and TCW, and was designed to 
highlight the un-vegetated spectral signatures and separate them from all other vegetation 
signatures. The DI metric is relatively insensitive to the changes in solar geometry or the 
bidirectional reflectance distribution function (BRDF) between scenes, and reduces the effect of 
phenologic variability among images obtained on different dates (Masek et al., 2008). Areas with 
disturbances have higher DI values. The advantage of the DI is that typically depends on the study 
region and only requires a threshold to be set. The DI has been used to successfully detect areas of 
forest disturbance in a wide range of forest biomes and DI=3.0 has been shown to provide the most 
accurate disturbance map (Baumann et al., 2014). The TCA introduced by Powell et al. (2010) 
condenses the information of the TCG/TCB ratio to a single value and fundamentally stands for the 
proportion of vegetated to non-vegetated areas; changes in the proportion of vegetated to non- 
vegetated areas result in concomitant changes in TCA values. Areas of dense vegetation have high 
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TCA values, while areas of bare soil with no vegetation have negative TCA values (Gómez et al., 
2011). The process indicator (PI) is a derivative of the TCA values, and represents the rate of change 
of TCA values over time. It is appropriate for detecting continuous subtle changes such as natural 
succession (Gómez et al., 2011). However, to the best of the authors’ knowledge, few studies had 
been conducted on sand dune encroachment using the TCT technique. 

The goal of this study is to characterize the desertification process using the TGSI with the three 
indices derived from the TCT. Another goal is to use remote sensing data from 2000 to 2015 to 
assess changes and trends in desertification in the MP as a whole based on the dune and sandy land 
classification of the decision tree (DT). In this study, thirteen sand dunes (deserts) and sandy lands 
were selected as the study sites to indicate the desertification processes of the MP as a whole. 


2 Materials and methods 


2.1 Study area 


The MP is located in the hinterland of temperate Asia and has a total population of about 2.8x107 
(0.3x107 for Mongolia and 2.5x107 for IMAR in 2015; Fang et al., 2015). It is dominated by steppe 
vegetation, including meadows, and typical and desert steppe vegetation (Zhao et al., 2015). It has 
a hyper-continental climate. The highest annual mean temperature appears in the desert steppe 
(1.2°C in Mongolia and 5.3°C in the IMAR) and the lowest annual mean temperature is —1.1°C in 
the desert steppe of Mongolia and 1.9°C in the desert steppe of IMAR (Fang et al., 2015; Zhao et 
al., 2015). The wind blows from between the northwest and the west for most of the year (Yao et 
al., 2007; Li, 2011; Karnieli et al., 2014; Lamchin et al., 2016). There are deserts of Badain Jaran, 
Tengger, Ulan Buh, and Hobgq, and sandy lands of Mu Us, Otindag, Horqin, and Hulun Buir in the 
IMAR (Zha and Gao, 1997; Kawamura and Akiyama, 2010). The Gobi is a vast, large desert region 
situated on the MP, and covers parts of northern IMAR and southern Mongolia. This study was 
carried out in Khar Nuur Sand Dune (Site 1), Ereen Nuur Sand Dune (Site 2), Tsagan Nuur Sand 
Dune (Site 3), Khongoryn Els Sand Dune (Site 4), Khar Khorin and Bayan Gobi (Site 5), Badain 
Jaran Desert (Site 6), Tengger Desert (Site 7), Ulan Buh Desert (Site 8), Hobq Desert (Site 9), and 
Mu Us Sandy Land (Site 10), Otindag Sandy Land (Site 11), Horqin Sandy Land (Site 12), and 
Hulun Buir Sandy Land (Site 13) in the MP (Fig. 1). 
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Fig. 1 MODIS mosaic image (2014) for Mongolia Plateau (MP) with locations of deserts and sandy lands. Sites 
1-13 are located in the Khar Nuur Sand Dune, Ereen Nuur Sand Dune, Tsagan Nuur Sand Dune, Khongoryn Els 
Sand Dune, Khar Khorin and Bayan Gobi, Badain Jaran Desert, Tengger Desert, Ulan Buh Desert, Hobq Desert, 
Mu Us Sandy Land, Otindag Sandy Land, Horgin Sandy Land, and Hulun Buir Sandy Land, respectively. 


2.2 Data and methodology 


MODIS Nadir BRDF adjusted reflectance (NBAR) data were obtained from the United States 
Geological Survey (USGS). Compared with raw reflectance, NBAR data minimizes artifacts in the 
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sample related to variable geometry, and has seven spectral bands comparable to those of Landsat 
TM, and is a logical choice for tasseled cap extension (Lobser and Cohen, 2007). For this study, 
the most recent collection 5 MODIS MCD43A4 16-day NBAR product at a spatial resolution of 
500 m acquired in August from 2000 to 2015 was freely downloaded from the global visualization 
interface (GloVis; http://glovis.usgs.gov/). Through mosaicking, projection transformation and 
cloud removal, a cloud-free dataset covering the study area was analyzed for the period 2000-2015. 
The TCB, TCG, and TCW are derived from MODIS MCD43A4 data using the tasseled cap 
coefficients given by Lobser and Cohen (2007). 

The TGSI was proposed based on field-based spectral reflectance measurements and laboratory- 
based physical analysis of topsoil grain composition, which is positively correlated (R?=0.7387) 
with fine sand content (Xiao et al., 2006). According to the results from Landsat TM/ETM+ (Xiao 
et al., 2006), a TGSI value close to zero corresponds to vegetation and water bodies. Sometimes, it 
can even be negative, and areas covered by fine sand (i.e., deserts) have TGSI values near 0.20. 
The TGSI can be estimated using Equation 1. 


(Roi — Bos) 
(Roi + Bos + Goa) ` 
where Ri, Goa, and Bss are the red, green, and blue bands of the MCD43A4 data, respectively. 


TGSI = (1) 


2.3 Classification of sand dunes and sandy lands 


The DT is applied to classify the sand dunes and sandy lands, which employs a flowchart-like 
process to perform multistage classification based on a series of splitting thresholds. DT 
classification has been used to assess the status of desertification (Xu et al., 2009; Lamchin et al., 
2016). With no ground truth data, statistical criteria combined with the opinion of experts were 
used to detect the sand dunes and sandy lands in the high resolution images acquired for the 
corresponding years from Google Earth using the variables TCA, DI, and TGSI (Fig. 2). In this 
study, the land desertification was classified into six grades: none (zero desertification), low sandy 
land (LSL, low desertification), coarse sandy land (CSL, medium desertification), fine sandy land 
(FSL, high desertification), coarse sand dune (CSD, severe desertification), and fine sand dune 
(FSD, severe desertification). After the DT classification, the masks (the rectangle buffers) were 
created to constrain calculations and analysis to the thirteen study sites. 


TCA, DI, TGSI 
0.0<TCA<0.1 TCAS0.0 


y Y 


DI>3.5 | | DI<3.5 | DI23.5 é| DI<3.5 
yooo y —— 


y 4 


Y 
| TGSIz0.26 || 0.21<TGSI<0.26 | Non —0.05<TCA<-0.02 | | Non 
pema amen 0.21<TGSI<0.26 TGSI>0.26 


Fine sandy DI>5.0 | DI<5.0 | 
land Y y 
Coarse Low sand: Coarse Fine sand 
sandy land land sand dune dune 


Fig. 2 Decision tree for classifying sandy lands using tasseled cap angle (TCA), disturbance index (DI), and 
topsoil grain size index (TGSI) 


DI>4.0 


2.4 Accuracy assessment 


A detailed accuracy assessment of DT classification could not be carried out due to insufficient 
field data. Informal assessment of accuracy was done by comparing the DT classifications to 
visualization interpretation of MCD43A4 images, some field photos, and Landsat images and high 
spatial resolution images acquired for periods corresponding to the classification imagery, and the 
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previous results (Hu et al., 2002; Xu et al., 2009; Guo et al., 2010; Li, 2011; Zhang et al., 2012; 
Duan et al., 2014; Lamchin et al., 2016; Wang et al., 2017). Because there were not enough field 
data as the ground truth verification values, 1,4593 random sample pixels (one percent pixels of 
the thirteen sites) were generated from the DT classification results using the random sampling 
method respectively. After the random sample pixels were overlaid on the same period Landsat 
TM/ETM+ and OLI images covered the thirteen sites, the correct and incorrect classification pixels 
were determined by visual interpretation with the high resolution images from the Google Earth. A 
confusion matrix was used for accuracy checking of non-desertification and desertification by 
produce accuracy, user accuracy, overall accuracy and Kappa coefficient. 


3 Results 


3.1 TCA, DI and TGSI 


Figure 3 shows the dynamic changes of TCA (Fig. 3a), DI (Fig. 3b) and TGSI (Fig. 3c) values of 
the thirteen sites over the last 16 years. At Site 1, the mean TCA values showed a decreasing trend 
during the periods 2000-2002, 2004-2008, and 2010-2012. The mean DI values showed an 
increasing trend during the periods 2004—2008 and 2010-2012. The mean values of TGSI of time 
series images from 2000 to 2015 were between 0.203 and 0.229. Compared with Site 1, extensive 
desertification occurred at Site 2. The mean TCA values showed a decreasing trend during the 
periods 2000-2002, 2005-2007, and 2010-2012. The mean DI values showed an increasing trend 
during the periods 2007—2009 and 2010-2012, indicating a fluctuating change in the conditions of 
land surface. The mean values of TGSI of time series images from 2000 to 2015 were between 
0.217 and 0.232, which showed that the dominant soil particle size of the topsoil had a fluctuating 
increasing trend of fine sand from 2000 to 2006, and a fluctuating decreasing trend of fine sand 
from 2006 to 2015. At Site 3, eleven years out of the last 16 years had negative average TCA values, 
indicating that a low proportion of vegetated to non-vegetated areas were present at Site 2 over the 
last 16 years. In particular, from 2004 to 2010, all the average TCA values were negative, indicating 
a persistent desertification. Compared with Site 1 and Site 2, a high proportion of the desertification 
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Fig. 3 Dynamic changes of TCA, DI and TGSI values of the thirteen sites from 2000-2015 
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occurred at Site 3. The mean DI values showed an increasing trend from 2003 and 2011 to 2005 
and to 2013, and a decreasing trend during the periods 2001—2003, 2005-2007, and 2013-2015, 
respectively. The mean values of TGSI of time series images from 2000 to 2015 were between 
0.213 and 0.225. At Site 4, the average TCA values were negative in the fifteen years out of the last 
16, indicating that a low proportion of vegetated to non-vegetated areas were presented at the site 
during these years. In particular, the average TCA values were close to —0.05 in 2001, 2002, 2005, 
2009, and 2010, indicating a lower proportion of vegetated to non-vegetated areas existed during 
these years. During the last 16 years, the mean DI values were near 1.0, indicating a low proportion 
of sand dunes and sandy lands present at Site 4. The mean values of TGSI of time series images 
from 2000 to 2015 showed a slightly increasing trend during the periods 2000-2003 and 2010- 
2012, and a slightly decreasing trend from 2003 to 2005. At Site 5, the average TCA values were 
positive at all time periods because sand dunes and sandy lands accounted for a small portion of 
the landscape at the site where the mean TCA values showed a decreasing trend during the periods 
2000-2002 and 2012-2014, and an increasing trend during the periods 2004—2006, 2007-2009, 
and 2010-2012. The average DI values werer negative in the twelve years out of the last 16 years, 
indicating a low proportion of sand dune and sandy lands at Site 5. The mean values of TGSI of 
time series images from 2000 to 2015 were 0.200—0.226, and showed a slightly increasing trend 
during the periods 2002-2004 and 2006-2008, and a slightly decreasing trend during the periods 
2004-2006 and 2010-2013. 

At Site 6 in the Badain Jaran Desert, the average TCA values were negative at the whole period 
2000-2015, indicating a high proportion of the presence of non-vegetated to vegetated areas. The 
average TCA values were all near —0.05 and the standard deviations of TCA values were relatively 
low for the entire period, indicating a relative stabilization of the landscapes. The mean DI values 
were highly positive for the whole period, indicating a relative high temperature, low moisture and 
sparse vegetation compared to the five abovementioned sites. The mean values of TGSI of time 
series images from 2000 to 2015 were 0.215—0.235. Unlike the five abovementioned sites, for most 
image dates, the distribution of pixels with different TGSI values was bimodal, indicating that the 
soil particles were sorted by size into the fine and coarse grains, with the majority of the particles 
being coarse grains. At Site 7 in the Tengger Desert, the average TCA values were negative for the 
period 2000-2015 except for the year 2012, indicating a high proportion of non-vegetated to 
vegetated areas. The mean DI values were highly positive for the period, indicating a relative high 
temperature, low moisture and sparse vegetation in the desert. The mean values of TGSI of time 
series images from 2000 to 2015 were 0.241—0.260. Compared with the six abovementioned sites, 
the mean TGSI values were the highest, indicating more fine sands in the desert. At Site 8, the 
mean values of TCA were less than 0.01, indicating a high proportion of non-vegetated to vegetated 
areas in the Ulan Buh Desert. The mean DI values were highly positive for the study period, 
indicating a relative high temperature, low moisture, and sparse vegetation in the desert. The mean 
values of TGSI of time series images from 2000 to 2015 were 0.217—0.237. It is interestingly noted 
that compared with the Tengger Desert, the mean values of TGSI were more like that of the Badain 
Jaran Desert. At Site 9, the mean values of TCA were more than 0.01 because sand dunes and sandy 
lands accounted for a small portion of the landscape. In the twelve years out of the last 16 years, 
the average DI had positive values, indicating a relatively high temperature and low moisture in the 
Hobg Desert. The mean values of TGSI of time series images from 2000 to 2015 were 0.207—0.220. 
The distribution of pixels with different TGSI values was bimodal, indicating that the soil particles 
were sorted by size into fine and coarse grains, with a majority of the grains being coarse grained. 
At Site 10, the mean TCA values showed an increasing trend from 2005 to 2008 and a decreasing 
trend during the periods 2008—2011 and 2012-2015, indicating a decreasing change with periodical 
fluctuations in the proportion of vegetation to non-vegetation. The mean DI values were negative 
at all time periods, indicating a relatively high vegetation cover in the Mu Us Sandy Land. The 
mean values of TGSI of time series images from 2000 to 2015 were 0.202—0.235. At Site 11, the 
average TCA values were more than 0.15 for the study period, indicating a relatively high 
vegetation cover in the Otindag Sandy Land. The mean TCA values showed an increasing trend 
during the periods 2001-2004 and 2009-2012, indicating a periodically fluctuating increasing 
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change in the proportion of vegetated to non-vegetated areas. The mean TCA values showed a 
decreasing trend from 2012 to 2014, indicating a slow decreasing change in the proportion of 
vegetated to non-vegetated areas. Fifteen of the last 16 years had low positive values of the average 
DI which may have been due to the presence of water bodies and sparse vegetation in the Otindag 
Sandy Land. The mean values of TGSI of time series images from 2000 to 2015 were between 
0.180 and 0.218, and showed an increasing trend during the periods 2004—2007, 2008-2010, and 
2011-2014, indicating a periodically fluctuating increase in the proportion of fine grained to coarse 
grained particles of the topsoil. At Site 12, the minimum mean TCA of time series images from 
2000 to 2015 was in 2009, the average TCA values were more than 0.18 for the whole period, 
indicating a relatively high vegetation cover in the Horqin Sandy Land. The mean TCA values 
showed an increasing trend from 2009 to 2013, indicating a periodically fluctuating increasing 
change in the proportion of vegetated to non-vegetated areas. The mean TCA values showed a 
decreasing trend during the periods 2001-2003, 2004—2007, and 2013-2015, indicating a slow 
decreasing change in the proportion of vegetation to non-vegetation. Seven years out of the last 16 
years had low negative values of the average DI, which may be caused by sparse water bodies and 
vegetation in the Horqin Sandy Land. The mean values of TGSI of time series images from 2000 
to 2015 were 0.192—0.213. At Site 13, the average TCA values were more than 0.12 for the study 
period, indicating a relatively high vegetation cover in the Hulun Buir Sandy Land. The mean DI 
values were negative for the period, which may be due to the presence of sparse water bodies and 
vegetation in the Hulun Buir Sandy Land. The mean values of TGSI of time series images from 
2000 to 2015 were 0.167—0.203. Compared with Horgin Sandy Land, the coarse sand content of 
the topsoil was more than that in the Hulun Buir Sandy Land. 

Although the mean TGSI of all thirteen sites showed a slight fluctuation, the mean TGSI values 
did not show a significant change over the period 2000-2015, indicating a relative stabilization in 
the soil grain size of the topsoil for the MP. 


3.2 Desertification change processes-PI 


The average PI value represents the global status of change, and can be used to discover trends in 
the study area through evaluation of the average PI at consecutive dates (Gómez et al., 2011). For 
most image dates, the distribution of pixels with different PI values was unimodal. The PI values 
were divided into five groups to explore the desertification change patterns. The stable group with 
near zero PI described areas with no changes in the proportion of non-vegetated to vegetated areas. 
According to the mean and standard deviation of PI at all dates, group thresholds (approximately 
the mean+two standard deviations) were determined: the PI value of the slow increase group (i.e., 
under slow desertification) ranged from —0.05 to 0.00, and the PI value of the slow decrease group 
(i.e., under slow restoration) ranged from 0.00 to 0.05. The pixels with PI values less than —0.05 
were sorted into the fast increase group (i.e., in a rapid desertification), and the pixels with PI 
greater than 0.05 were sorted into the fast decrease group (i.e., in a rapid restoration). 

Figure 4 shows the dynamic changes of the PI values of the thirteen sites over the last 16 years. 
At Site 1, nine years out of the last 16 years had low negative average PI values. In particular, the 
average PI values were negative from 2004 to 2008, indicating a slow increasing rate in the 
proportion of non-vegetated to vegetated areas, i.e., desertification. The mean PI values in 2004, 
2005, 2007, 2008 and 2012 were close to zero, indicating a relative stabilization during these years. 
At Site 2, ten years out of the last 16 years had negative average PI values. The average PI values 
were negative from 2004 to 2008, indicating a slow increasing rate in the proportion of non- 
vegetated to vegetated areas, i.e., desertification. It should be noted that a significant change in 
restoration occurred in 2002 and 2010. The mean PI values in 2005, 2007, 2008, 2012, 2013 and 
2014 were close to zero, indicating relative stability during these years. At Site 3, eight years out 
of the last 16 years had negative values in the average PI. The average PI values were positive from 
2009 to 2011, indicating a slow decreasing rate in the proportion of non-vegetated to vegetated 
areas, whereas the average PI values were negative from 2012 to 2014, indicating a slow increasing 
rate in the proportion of non-vegetated to vegetated areas. The mean PI values in 2003, 2005, 2008, 
2009, 2012, 2013 and 2014 were close to zero, indicating a relative stabilization during these years. 
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At Site 4, nine years out of the last 16 years had negative average PI values. The average PI values 
were negative from 2007 to 2009, indicating a slow increasing rate in the proportion of non- 
vegetated to vegetated areas. The mean PI values in 2003, 2005, 2007, 2012, 2013 and 2014 were 
close to zero, indicating relative stability during these years. At Site 5, ten years out of the last 16 
years had positive average PI values. The average PI values were positive from 2008 to 2012, 
indicating a slow decreasing rate in the proportion of non-vegetated to vegetated areas. The mean 
PI values in 2007, 2009, and 2010 were close to zero, indicating relative stability of the landscape 
during these years. The mean PI values in 2001 and 2013 were near —0.05, indicating a relative 
rapid desertification during these years. The mean PI values in 2005 and 2011 were near 0.05, 
indicating a relative rapid change toward restoration during these years. For sites 1—5, the pixels in 
the fast increase and decrease groups were relatively infrequent during the study period, and the 
pixels in the slow increase and decrease groups were normally frequent. 
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Fig. 4 Dynamic changes of the process indicator (PI) values of the thirteen sites from 2000-2015 


At Site 6, nine years out of the last 16 years had positive average PI values, indicating a slow 
decreasing rate in the proportion of non-vegetated to vegetated areas. However, the mean PI values 
were close to zero, indicating a relative stability of the Badain Jaran Desert over the last 16 years. 
At Site 7, nine years out of the last 16 years had negative average PI values, indicating a slow 
increasing rate in the proportion of non-vegetated to vegetated areas. However, the mean PI values 
were all close to zero, indicating relative stability of the Badain Jaran Desert over the last 16 years. 
At Site 8, ten years out of the last 16 years had negative average PI values, indicating a slow 
increasing rate in the proportion of non-vegetation to vegetation. However, the mean PI values were 
all close to zero, indicating relative stability of the Ulan Buh Desert during the last 16 years. At 
Site 9, eight years out of the last 16 years had negative values in the average PI, indicating a slow 
increasing rate in the proportion of non-vegetated to vegetated areas. The mean PI values in 2003, 
2005 and 2008 were close to zero, indicating relative stability during these years. At Site 10, nine 
years out of the last 16 years had the negative average PI values, indicating a slow increasing rate 
in the proportion of non-vegetated to vegetated areas. The mean PI values were close to zero during 
the periods 2003-2005, and 2007-2010, indicating a relative stabilization during these years. At 
Site 11, nine years out of the last 16 years had the negative average PI values, indicating a slow 
increasing rate in the proportion of non-vegetated to vegetated areas. The mean PI values were 
close to zero in 2001, and 2005 to 2007, indicating a relative stabilization during these years. At 
Site 12, nine years out of the last 16 years had positive values in the average PI, indicating a slow 
decreasing rate in the proportion of non-vegetated to vegetated areas. At Site 13, ten years out of 
the last 16 years had positive average PI values, indicating a slow decreasing rate in the proportion 
of non-vegetated to vegetated areas. The mean PI values were close to zero in 2006, 2007, 2009, 
2011, and 2012, indicating relative stability during these years. For sites 6—13, the pixels in the 
slow increase and decrease group were more common during the majority of the study period. 
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3.3 Assessment of desertification 


The 500 m pixels of MCD43A4 products represent a mixture of land cover classes that influence 
the spectral ranges, and other objects (such as the Gobi desert, and barren lands) with similar TCA, 
DI and TGSI as sand dunes and sandy lands, which could also create some errors. Through the 
buffer set on the study site, the DT classification errors could be reduced partly. The same period 
Landsat TM/ETM-+ (fifty-two Landsat 5 TM and five Landsat 7 ETM+ images) and OLI images 
(nineteen Landsat 8 OLI images) covered the thirteen sites (Site 1, path/row 139/26 and 138/27; 
Site 2, path/row 139/27 and 138/27; Site 3, path/row 134/29; Site 4, path/row 133/30; Site 5, 
path/row 133/27; Site 6, path/row 133/32 and 132/32; Site 7, path/row 130/33; Site 8, path/row 
130/32; Site 9, path/row 129/32 and 128/32; Site 10, path/row 128/32 and 128/33; Site 11, path/row 
124/30 and 125/30; Site 12, path/row 121/30 and 122/30; Site 13, path/row 124/26), which were 
used to evaluate the DT classification. As Table 1 showed, the overall Kappa coefficient and overall 
accuracy values of the DT classification results were 0.7759 and 91.71%, 0.7659 and 90.71%, 
0.7737 and 91.71%, 0.7649 and 90.44% in 2000, 2005, 2010, and 2015, respectively, which 
reflected the overall classification situation. TGSI is a key factor for classifying the desertification 
into CSL, FSL, CSD, and FSD. Compared with the grain size analysis results of the deposits in the 
Badian Jaran-Tengger Deserts (Li, 2011), it showed that TGSI was a good indicator for the grain 
size analysis of the surface deposits, which indirectly indicated that the desertification classification 
of this paper was credible. In addition, the desertification maps derived from the DT classification 
were comparable to the desertification map of Mu Us Sandy Land in 2000 and 2010 (Xu et al., 
2009; Wang et al., 2017), Horqin Sandy Land in 1999, 2000, and 2010 (Hu et al., 2002; Zhang et 
al., 2012; Duan et al., 2014), Hulun Buir Sandy Land in 2000 and 2006 (Guo et al., 2010), and Khar 
Khorin and Bayan Gobi in 2002 and 2011 (Lamchin et al., 2016), which could provide the reliability 
of the classes of interest and could be employed to assess the desertification on a regional scale. 
Due to the excessive number of tables and figures for the thirteen sites for the last 16 years, we will 
only show the desertification map of the thirteen sites in 2015 (Fig. 5), and for the annual changes 
of the main sand dune of Site 3 from 2000 to 2015 (Fig. 6). 


Table 1 Accuracy analysis for desertification for the thirteen study sites in 2000, 2005, 2010, and 2015 


Number of pixels Number of pixels 

Classes Year Year 

Sand dune Non Total Sand dune Non Total 
Sand dune 3445 592 4037 2920 534 3454 
Non 803 9753 10,556 676 10,463 11,139 
Total 4248 10,345 14,593 3596 10,997 14,593 
Producer accuracy 2015 81.10% 94.28% 2010 81.20% 95.14% 
User accuracy 85.34% 92.39% 84.54% 93.93% 
Overall Kappa coefficient 0.7649 0.7737 
Overall accuracy 90.44% 91.71% 
Sand dune 3306 624 3930 3398 770 4168 
Non 732 9931 10,663 543 9882 10,425 
Total 4038 10,555 14,593 3941 10,652 14,593 
Producer accuracy 2005 81.87% 94.09% 2000 86.22% 92.77% 
User accuracy 84.12% 93.14% 81.53% 94.79% 
Overall Kappa coefficient 0.7659 0.7759 
Overall accuracy 90.71% 91.00% 


Note: Sand dune includes low sandy land (LSL, low desertification), coarse sandy land (CSL, medium desertification), fine sandy land 
(FSL, high desertification), coarse sand dune (CSD, severe desertification) and fine sand dune (FSD, severe desertification); Non, no 
sandy land. 


From 2000 to 2015, it is clear that most of the study sites showed fluctuating trends of 
desertification, and not a monotonic increase. The area of desertification of Site 1 (20.6% of the 
total area), Site 2 (34.0% of the total area), Site 3 (51.1% of the total area), Site 6 (65.7% of the 
total area) and Site 7 (68.1% of the total area) was the highest in 2012. The area of desertification 
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of the Site 4 (18.3% of the total area), Site 5 (5.6% of the total area), Site 8 (55.2% of the total 
area), Site 9 (15.6% of the total area), Site 10 (13.4% of the total area) , Site 11 (21.6% of the total 
area), Site 12 (20.0% of the total area), and Site 13 (3.8% of the total area) was at the highest in 
2001, 2002, 2003, 2003, 2000, 2005, 2009, and 2004, respectively. 

Figure 6 shows the annual changes of the main sand dune in Site 3 over the last 16 years. In 
2001, total area of desertification was at its maximum, indicating the occurrence of a severe drought 
in 2001. The total desertification area and the area of coarse sand dune showed an increasing trend 
from 2003 to 2005, and decreasing trends in 2001—2003, 2005-2007, 2008-2011, and 2012-2015; 
and the overall trend was slow down from 2000 to 2015 (Fig. 6a). Figure 6b indicated a fluctuating 
increasing trend of the area of low sandy land, fine sandy land and fine sand dune. The southwest 
and northeast parts of the main sand dune underwent fluctuating changes over the last 16 years. 
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Fig. 6 Annual area changes of the main sand dune in Site 3 for the period 2000-2015. Total, total desertification 
area; CSD, coarse sand dune (severe desertification); LSL, low sandy land (low desertification); FSD, fine sand 
dune (severe desertification); FSL, fine sandy land (high desertification). 


For sand dunes and deserts, the coarse sand dunes and the fine sand dunes occupied the majority 
of the desertified area. The areas of the coarse sand dune and the fine sand dune were at their highest 
in 2008 (6.2% of the total area) and 2013 (1.5% of the total area) at Site 1, in 2009 (20.4% of the 
total area) and 2013 (3.3% of the total area) at Site 2, in 2008 (42.4% of the total area) and 2005 
(0.9% of the total area) at Site 3, in 2001 (15.4% of the total area and 2.6% of the total area 
respectively) at Site 4, in 2003 (54.2% of the total area) and 2012 (8.7% of the total area) at Site 6, 
in 2007 (41.9% of the total area) and 2006 (26.1% of the total area) at Site 7, in 2003 (51.3% of the 
total area) and 2013 (12.8% of the total area) at Site 8, in 2001 (8.2% of the total area) and 2003 
(2.2% of the total area) at Site 9, respectively. For the sandy lands, the low sandy land, coarse sandy 
land, and fine sandy land occupied the majority of the desertified area. The areas of the low sandy 
land and the fine sandy land were the highest in 2002 (5.3% of the total area) and 2004 (1.7% of 
the total area) at Site 5, in 2000 (7.7% of the total area and 3.8% of the total area) at Site 10, in 
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2014 (7.8% of the total area) and 2004 (1.7% of the total area) at Site 11, respectively. However, 
the areas of the low sandy land and the coarse sandy land were the highest in 2009 (14.6% of the 
total area) and 2000 (5.7% of the total area) at Site 12, and in 2004 (3.4% of the total area and 0.3% 
of the total area at Site 13, respectively). 


4 Discussion 


The TCA value describes the proportion of vegetated to non-vegetated areas. The TCA used in this 
study is a valuable parameter for the discrimination of sand dune, sandy land, and vegetation. In 
general, the TCA values of sand dunes were less than zero while the TCA values of the sandy lands 
were 0.0-0.1, indicating a better vegetation cover. The minimum mean TCA values of sites 1, 5, 
and 6 occurred in 2002; the minimum mean TCA values of sites 2, 11, and 12 occurred in 2009; 
and the minimum mean TCA values of sites 3, 4, 7, 8, 9, 10, and 13 occurred in 2001, indicating 
the occurrence of a severe drought and lower soil moisture content during these years (Dorjsuren 
et al., 2016; Lamchin et al., 2016). It is worth examining why the minimum mean TCA values of 
the different sites occurred in different years, which has thus far been attributed to different climate 
change rates between the different sites. Compared with the TCA, the DI is created using TCB, 
TCG, and TCW, which reduced the effects of moisture and vegetation simultaneously. In some 
cases, the DI is better able to identify sand dunes from the background than the TCA. The TGSI 
reports the grain size of the surface in a defined area. In this study, the TGSI values of more than 
0.21 were classified into the sand dune or sandy land categories, unlike the studies using Landsat 
images in the Hogno Khaan protected area in Mongolia and the arid to semi-arid regions of Inner 
Mongolia (Xiao et al., 2006; Lamchin et al., 2016). This difference is likely due to the use of 
different sensors, although the spectral bands of the MCD43A4 have been adjusted to be 
comparable to those of Landsat TM. Compared with the DI, TGSI can better describe the structures 
of the sand dunes and sandy land, such as the coarse and fine sand dune, and the coarse and fine 
sandy land, which is important for studying the sand dune movement, sand shift, and so on. 
However, in some cases, it cannot adequately differentiate the sandy land from the background. 
Compared with the other sand dunes of the study area, the Tengger Desert had higher DI and TGSI 
values, indicating a brighter and finer surface. Therefore, the integration of the advantages of TCA, 
DI, and TGSI could better identify and assess the desertification. 

The PI is a valuable parameter used in assessing trends and processes of desertification change 
that might be missed by an exclusively bi-temporal change detection approach. It should be noted 
that the maximum mean PI did not correspond to the year with the maximum precipitation, which 
indicated that the restoration might have reached its environmental upper limit before the year with 
the maximum precipitation. During the last 16 years, the mean PI values of Badain Jaran, Tengger, 
and Ulan Buh deserts were all close to zero, indicating a relative stabilization of these deserts. The 
mean PI values were negative at sites 1, 2, 3, 4, 9, 10, and 11 in more than half of the study period, 
indicating a slow increase in desertification, whereas the mean PI values were positive at sites 5, 
12, and 13, indicating a slow decrease in desertification. It has not been determined whether the 
desertification is permanent or whether it is only an expression of seasonal variability (Sternberg 
et al., 2011; Eckert et al., 2015). However, our study identified parallel processes of ongoing climate 
variability and fluctuations in desertification, indicating a three or four-year periodically fluctuating 
change rather than a continuous increase or decrease in the desertification, which may be reflected 
by changes in precipitation, soil moisture, and vegetation restoration projects in the MP (Sternberg 
et al., 2011, 2015; Dorjsuren et al., 2016; Huang, 2017; Yu et al., 2017). The Horqin Sandy Land 
in particular, showed random alternations between a slow fluctuated increase and decrease for all 
the time periods, indicating the impacts from the local ecological restoration projects (Zhang et al., 
2012). 

For the Khar Nuur, Ereen Nuur, Tsagan Nuur, and Khongoryn Els sand dunes, and the Badain 
Jaran, Tengger, Ulan Buh and Hobq deserts, the fine sand dunes were mainly distributed in the east 
and southeast parts of the sites, indicating that the fine sand shifts with the predominant winds. Our 
study showed that compared with the other sand dunes and deserts, the fine sand dunes occupied 
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the majority of the Tengger Desert, as is evidenced by the results of studies based on the grain size 
of deposits in the Badain Jarden-Tengger Desert (Li, 2011), which indirectly indicate that the DT 
classification of this study was receivable, and the coarse sandy land occupied the majority of the 
Horqin Sandy Land. Our study showed that sandy lands around the dunes or deserts had finer sand 
particles than the coarse sandy land, for example, comparisons of sites 1 and 12 show the 
differences in topsoil grain sizes between the west and the east of the MP. Further studies need to 
be conducted in the future. 


5 Conclusions 


Desertification is a complex process driven by climate change and anthropogenic impacts. In this 
study, we attempted to monitor the desertification changes in the MP from 2000 to 2015 through 
analysis of the thirteen sites, including sand dunes and sandy lands, using MODIS MCD43A4 data. 
Although with no ground truth data were available to carry out a detailed accuracy assessment of 
the DT classification, visual interpretation, and high resolution images from Google Earth acquired 
for the same periods, together with previous research results showed that the DT classification 
could be employed to assess the desertification at the regional scale. 

It demonstrated that the TCA is a valuable parameter for the assessment of sand dunes, sandy 
lands, and vegetation. In some cases of this study, the DI was better able to identify sand dunes 
from the background than the TCA. The TGSI was able to better describe the structures of the sand 
dune and sandy land, such as the coarse and fine sand dunes, and the coarse and fine sandy lands, 
which were important for studying sand dune movements, and sand shift. The PI provided a view 
of the trends and processes of desertification change that might be missed with only a bi-temporal 
change detection approach. The integration of the advantages of the TCA, DI, TGSI, and PI could 
better identify and assess the desertification. During the last 16 years, the deserts of Badain Jaran, 
Tengger, and Ulan Buh were relatively stable. The Khar Nuur, Ereen Nuur, Tsagan Nuur, 
Khongoryn Els, and Hobq deserts, and the Mu Us, and Otindag sandy lands showed a gradual 
increase in desertification, whereas the Bayan Gobi, Horqin, and Hulun Buir sandy lands showed 
a gradual decrease in desertification. 

At the moment, various limitations remain, including the mixed pixel effects and threshold 
selection. Although evaluating relative changes in TCA, DI, TGSI, and PI did not require 
calibration, and MODIS MCD43A4 had standardized reflectance to a nadir view, normalization 
between the TCA, DI, TGSI, and PI values should be studied in the future. Further studies to link 
the TCA, DI, TGSI, and PI values with the desertification characteristics are recommended to set 
the thresholds and improve the assessment accuracy. 
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